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Abstract. The classical limit problem of quantum mechanics is revisited on 
the basis of a scheme that enables a quantitative study of the way the quantum- 
classical agreement emerges while going through the intermediate mass range 
between the microscopic and the macroscopic domains. As a specific application 
of such a scheme, we investigate the classical limit of a quantum time distribution 
- an area of study that has remained largely unexplored. For this purpose, we 
focus on the arrival time distribution in order to examine the way the observable 
results pertaining to the quantum arrival time distribution which is defined in 
terms of the probability current density gradually approach the relevant classical 
statistical results for an ensemble that corresponds to a Gaussian wave packet 
evolving in a linear potential. 
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1. Introduction and the underlying basic scheme 

Over the years, the analysis of various aspects of the classical/macroscopic Umit of 
quantum mechanics has attracted considerable attention[l-10]. Broadly speaking, 
there are two distinct strands of investigations related to the classical limit of quantum 
mechanics. One direction of study has been to delineate the way 'classical-like 
behavior' can be obtained for any quantum mechanical micro-system under suitable 
conditions. The other line of study seeks to probe the macroscopic range of validity of 
quantum mechanics by examining as to what extent the quantum mechanical results 
in a suitably defined macroscopic regime agree with the corresponding results derived 
from classical mechanics. It is from the latter perspective that, in this paper, we 
investigate the quantum-classical correspondence for an observable time distribution 
- an issue that has been hitherto neglected in the context of the classical limit aspect 
of quantum mechanics. 

A particular aspect of the present paper is that the quantum-classical transition 
is probed in a quantitative way. For this purpose, in order to characterize the relevant 
macroscopic domain, we use 'mass' as the parameter that is varied to study the way the 
convergence of classical and the corresponding quantum results occurs by approaching 
the large mass values while going through the intermediate range. At this stage, some 
relevant remarks would be appropriate about the legitimacy of treating 'mass' as a 
parameter, instead of taking it to be an operator. 

First, note that the scope of our analysis is restricted within the nonrelativistic 
domain being based on the Galilean invariant Schrodinger equation for the spin-0 
particles. Now, given the Galilean invariance of the Schrodinger equation, one may 
recall an interesting theorem due to Bargmann[Tl] which states that, in nonrelativistic 
quantum mechanics, one cannot have a coherent superposition of states of different 
masses(for an elegant proof of this theorem, using a suitable sequence of Galilean 
constant velocity transformations, see, for example, Kaempffer[T2]). An insight into 
the physical justification for this theorem has been provided by Greenberger[13j 
who showed that this restriction arises essentially because that in the nonrelativistic 
domain, the 'coordinate time' does not differ from the 'proper time' (measured in the 
moving frame). One is, therefore, entitled to treat 'mass' as a parameter, as long as 
the study is restricted within the framework of nonrelativistic quantum mechanics. 

Next, coming to the conceptual basis of the quantum-classical comparison scheme 
that will be specifically used in this paper, we first note the following. While 
the predictions of the quantum mechanical formalism are verifiable pertaining to 
essentially an ensemble of particles [131 [H], classical mechanics can describe the 
properties of an ensemble of particles as well as of a single particle. Thus, the 
comparison between these two mechanics is operationally unambiguous provided one 
compares their statistical predictions for the dynamical evolutions of the same given 
initial ensemble. 

It is in this spirit that we adopt the scheme used in this paper [16] where the 
quantum and the classical evolutions are compared by starting from the same initial 
ensemble that has the specified position and momentum distributions obtained from 
a given wave function. While the quantum evolution is in accordance with the 
Schrodinger equation, the classical evolution of the given initial ensemble is calculated 
in terms of the classical phase space dynamics based on Liouvillc's equation. However, 
a critical point in the classical calculation is the following. The initial phase space 
distribution for an ensemble is not uniquely fixed even if the position and the 
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momentum distributions arc specified. But, since in classical mechanics, the time 
evolution of all the usual observable properties of an ensemble are determined by the 
initial positions and momenta which are mutually independent variables, an initial 
phase space distribution Do{xo,pQ,t — 0) evolving under classical dynamics can be 
written, in the simplest possible choice, as a product of the position and momentum 
distributions pertaining to a given initial wave function ^(xo,t ~ 0), given by 

Do{xo,Po,t = 0) = |*(xo,0)p |$bo,0)P (1) 

where the variables xq and po are the initial positions and momenta of the particles, 
and $(poi 0) is the Fourier transform of ^'(a:o: 0)- 

Based on this specific quantum-classical comparison scheme, the plan of this paper 
is as follows. In Section II we discuss the basic features of both the quantum and the 
classical procedures for defining the arrival time distribution using the probability 
current density. Here we may note that in recent years, the quantum mechanical 
distributions of various types of time like the tunneling time, arrival time, transit 
time, decay time, and so on have been widely studied; for useful reviews, see, for 
example, Muga et a/. [17] and Olkhovski et aZ.[T8]. In the light of this flourishing 
line of works, the classical limit aspect of such quantum time distributions deserves 
to be a germane area of study. To this end, in this paper, we initiate such an 
investigation by restricting our attention to the classical limit of a particular form of 
quantum arrival time distribution that is defined in terms of the probability current 
densitv [T7l \W[ \n\ - our analysis being contingent upon a specific scheme for the 
quantum-classical comparison, and is couched in terms of a Gaussian wave packet 
propagating in a linear potential, while such a study, in principle, can be extended for 
other forms of time-distributions, using wave functions of various types, and in the 
context of any other potential. 

In Section III, the quantum-classical correspondence of an arrival time 
distribution is treated in detail in terms of a general Gaussian wave packet (that 
does not correspond to the minimum value of the uncertainty product AxAp) which 
evolves in the presence of a one dimensional linear potential. The salient feature of 
this work is a quantitative study that is aimed at delineating the mass range over 
which the quantum results pertaining to the time distribution under consideration 
gradually concur with their classical counterpart - the representative relevant results 
for the mean arrival time and the associated fluctuation being given in section IV. In 
the concluding Section V some future directions of work are indicated. 

But, before proceeding further, for the sake of completeness, some remarks are 
in order to stress the conceptual inadequacy of the usual textbook definition of the 
classical limit of quantum mechanics given in terms of h 0. First, the notion 
that ?i ^ is 'small' has no absolute meaning because its value depends on the 
system of units [22]. Further, wave functions are, in general, highly nonanalytic in the 
neighborhood of the limit point h Q [23] . This results in the essential singularity of 
the quantum mechanically computed quantities at ?i — > . It is, thus, not possible to 
regard quantum mechanics as a perturbative extension of classical mechanics in the 
same sense as special relativity can be viewed as related to Newtonian mechanics by 
a convergent perturbation expansion in w/c[2j. Hence, the only sensible operational 
formulation of the h —^ classical limit condition would be to consider a dimensionless 
parameter of the form H/S << 1 where S is the 'action quantity' relevant to a given 
situation. But, then, within this approach, an clement of arbitrariness comes into play 
in the choice of the appropriate 'action quantity' to be used in any given example for 
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probing the classical limit of quantum mechanics. In contrast, the procedure adopted 
in our paper for studying the macrolimit of quantum mechanics by varying 'mass' as 
the relevant parameter is devoid of any such arbitrariness. 



2. Arrival time distributions in quantum and classical dynamics 

First, let us consider the quantum mechanical case. For simplicity, throughout this 
paper, we restrict the treatment to one spatial dimension. We begin with the non- 
relativistic quantum mechanical description of the flow of probability, expressed in 
terms of the position space distribution, that is governed by the continuity equation 
(derived from the Schrodinger equation) given by 

||*(x,i)p + V.J(x,t) = (2) 

The quantity J(x, t)=^(^'V** - **V*), called the probability current density, 
characterises this flow of probability. It is this current density that has been used in a 
number of studies to define the arrival time distribution for free particles [171 [TOl 
By interpreting the equation of continuity in terms of the flow of physical probability, 
in conjunction with using the Born interpretation for the squared modulus of the 
wave function as denoting the probability density, it has been suggested that the 
mean quantum arrival time of the particles reaching a detector located at x = X may 
be written as 

_ nj(x=x,t)M 

J-\Ji. = X,t)\dt 

whence the corresponding fluctuation (At)Q is given by the root mean square deviation 

iAt)Q = ^{tl)-^{tQ)\ 

The definition of the mean arrival time specified by Eq.Q is, however, not 
a uniquely derivable result within standard quantum mechanics. In fact, different 
schemes for defining the quantum arrival time distribution have been discussed in 
the literature; for example, using Kijowski's axiomatic approach[25], or by invoking 
the time-of-arrival operator method in conjunction with the POVM approaches], by 
constructing self-adjoint variants of the time-of-arrival operator [SJ, or by using the 
Bohmian causal model |28]. However, the ambit of the present paper is confined to 
only the probabihty current density approach [Tf t fT9 l [20] . Here it may be noted that 
in certain situations, the quantity J{X,t) can be negative during some time interval, 
even if the initial wave function has the positive momentum support - this is called 
the backflow effect [5^ . It is in order to take this effect into account that the modulus 
of the quantity J(X, i) (suitably normalised) is taken for specifying (ig) as given by 
Eq.Q. 

Next, we note that the Schrodinger probability current defined in terms of the 
continuity equation has an inherent ambiguity. This is because the continuity equation 
remains satisfied with the addition of any divergence free term to the probability 
current. On this point, Holland |30| has shown the uniqueness of the probability 
current for the spin-1/2 particles using the Dirac equation. On the other hand, for 
the spin-0 particles, using the Kemmer equation [31j . it has been demonstrated [32] 
that the non-relativistic limit of the Kemmer probability current is unique, whose 
expression turns out to be that of the Schrodinger probability current. Hence, for the 
spin-0 particles, the Schrodinger probability current can be used for computing the 
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arrival time distribution. Thus, even though the Schrddinger probabihty current is 
not directly observable, having no correspondence with an appropriate self-adjoint 
operator [171 I33| . it can have an observable manifestation for the spin-0 particles 
through the arrival time distribution. The latter is, in practice, a measurable quantity 
- this point being underscored in various experimental contexts involving the time- 
of-flight measurements concerning, for example, cold trapped atoms, with the 
quantum probability current being invoked in the relevant theoretical analysis [35]. 
Besides, several theoretical models of 'quantum clock' [211 [36] have been studied that 
bring out the empirical relevance of time distributions such as the arrival/transit time. 

Now, let us examine the classical procedure for computing the arrival time 
distribution. For this, a classical statistical ensemble of particles is described by 
the phase space density function D(x,p,t) . Consequently, the classical position 
and momentum distribution functions are respectively pc{x,t) = J D{x,p,t)dp and 
PciPi ^) = / D{x,p, t)dx, while D{x,p, t) satisfies the classical Liouville equation given 

by 

dD{x,p,t) ^ ^ dD{x,p,t) ^ ^ dD{x,p,t) 
dt dx dp 

Integrating the above equation with respect to p one gets 



dpc{x,t) ^ d_ 
dt dx 



—p{x,t)pc{x,t) 
m 



= (5) 



where p = J pD{x,p,t)dp/ J D{x,p,t)dp is the ensemble average of the momentum 
values of the individual particles. 

Defining v{x,t) = p{x,t)/m as the ensemble average of the individual velocity 
values, we get 

dpc{x,t) d 

-^— + -Jaix,t) = (6) 

where Jc{x,t) — pc{x,t)v{x,t). Thus. Eq.® can be regarded as the equation 
of continuity characterising the classical time evolution of a statistical ensemble of 
particles. Using the expressions for pc{x,t) and v{x,t)^ the classical probability 
current density can then be written as 

Jc{x,t) = — [ pD{x,p,t)dp (7) 



Given this statistical description, the mean classical arrival time is given by 

|Jc(x = X,t)\tdt 
^'^'^ ^ J,-\Jc{^ = X,t)\dt 
whence the corresponding fiuctuation {A.t)c is given by the root mean square deviation 



3. Quantum-classical correspondence for a 

non-minimum-uncertainty-product wave packet propagating in a linear 
potential 

In this section, we compare the quantum and the classical results for the position, 
momentum and time distributions by considering a general non-minimum-unccrtainty- 
product Gaussian wave packet propagating in a linear potential {V — Kx). 
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Here we take the initial wave function \E'(a;,0) and its Fourier transform 0) 
to be given by 

^2 



*(a;,0) 



{2nao^y/''VTTiC 



=exp 



Aao^{l + iC) 



+ ikx 



exp 



(l + ^C) 



(9) 
(10) 



where the group velocity of the wave packet u = Tik/m = p/m. 

Note that we have taken an initial Gaussian wave function ^'(a;,0) which is not 
a minimum uncertainty state, i.e., AxAp = (?i/2)vT+"C^ > Ti/2, where C is any 
real number - such a non-minimum-uncertainty-product state corresponds to what is 
known as a squeezed state |37] . In the presence of a linear potential, for such an initial 
wave function, the Schrodinger time evolved wave function ^(a;,<), and consequently 
the probability current density jQ{x,t) are respectively given by 



^{x,t) 



(2wo')i/yi+*(c- 



ht 



exp 



2m(To 



X exp ■ 



[x 



ut^\ief 



l + i(C- 



ht \ 



Kt 



ut , 



{u-—){x--) 



2j-3- 



iKH 
6mh 



(11) 



Kt HC 



Jqix^t) ^ PQ{x,t) in h 



^^){x~ut+\i^t') 



(12) 



where pQ{x,t) is the quantum mechanical position probability distribution function 
given by 



PQix,t) = |*(x,t)|^ 



1 



exp ■ 



2al{t) 



(13) 



where (7Q{t) = (7oyl + {C+ 2mao'^ )^ ^^"^ width of a quantum wave packet that 
corresponds to the position probability distribution at any given instant t. 

Next, we focus on calculating the probability current density J(.t, t) using the 
classical statistical evolution. For this purpose, in accordance with Eq.([T]), it is 
crucial that the initial phase space distribution Do(xo,pQ,t = 0) to be used for the 
classical calculations is fixed by the initial position and momentum distributions of 
the ensemble that are taken to be the same as the corresponding initial quantum 
distributions obtained from Eq.® and (jlOp respectively. Accordingly, the expression 
for i:»o(a;o,Po,* = 0) is given by 

Do{xo,po,0) = \^{xo,0)\^ |$(po,0)P (14) 



1 



exp 



2ao^ipo -p? 



[ 2ao2(l + C2) 

Now, in order to obtain the time evolved classical phase space density function 
D{x,p,t), we consider the classical Hamiltonian for the freely moving particles 



H = p^ /2m + A'.T, and Hamilton's equations given by a; = pot/m — ^^f^ 



xo 
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and p ^ Po — Kt. Then one can write Xq = x — p^t/m + \ ^t^ and Po ~ p + Kt. 
By substituting these values of xq and pq in the expression for DQ{xQ,p{i,Q) given 
by Eq. (fTi|) . we obtain the time evolved classical phase space distribution function 
D(x,p, t) given by 



2 



Now, substituting in Eq.Qthe expression for the time evolved phase space 
distribution function D{x,p, t) from Eg. ljlSp . the probability current density pertaining 
to this classical ensemble is given by 

( Kt {x^ut+\^f) 



1 ra ' 2mCTo^ 

in ' 2ma^{t) 



Jc{x, t) = pc{x, t) {u-^-^ + , ) (16) 



whence the position probability distribution for this classical ensemble is given by 

r if ix - ut + -—t^)'^\ 
pc{x,t)= D{x,p,t)dp= ^===xexp{-^ ^ 2 } (17) 



where ac (t) = ao^^J 1 + C"^ + is the time- varying width of the classical position 

distribution function at an instant t. 

Note that this spreading of the statistical distribution embodied in Eq. fTT)) 
ensues from the classical Liouville evolution, and is, in general, different from 
the corresponding quantum spreading of a wave packet((TQ(i) 7^ adt)), unless 
C ~ 0. This means that the quantum and the classical spreadings of the position 
distributions agree only if the initial Gaussian statistical distribution corresponds to 
the minimum- uncertainty-product, i.e., if initially, AxAp ~ Ti/2. It, therefore, needs 
stressing that such a spreading is not essentially a quantum mechanical property of 
a propagating wave packet, but is a generic feature associated with a time- varying 
position probability distribution, quantum or classical. While for the positive values 
of C, cTQ{t) is larger than crcit) for all times, for the negative values of C, <JQ{t) is 
always smaller that (Jc{t)- 

On the other hand, Egs. fT^ and clearly show that the quantum and the 
classical probability currents are, in general, not the same , i.e., jQ{x,t) ^ Jc{x,t). 
But, if one imposes the minimum-uncertainty-product condition C = 0, the quantum 
and the classical probability currents become the same, i.e., Jq{x, t) = Jc{x, t). Also, 
interestingly, this condition C = ensures an exact agreement sans any limiting 
condition between the quantum and the classical position probability distributions 
given by Eqs.(13) and (17) respectively. 



4. Results of some relevant quantitative studies and their implications 

In order to make a systematic study of the way an agreement emerges between the 
quantum and the classical probability currents given by Eqs.(12) and (16) respectively, 
thereby leading to a matching of the corresponding mean arrival times and their 
fluctuations for the non-minimum-uncertainty-product Gaussian wave function(C ^ 0) 
under consideration, we proceed as follows. First, for a fixed value of C, we compare 
the plots of the quantum and the classical probability currents by varying the values 
of the masses. As representative studies, we take C = 50 and the choices of the 
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Figure 1. The time variations of the quantum and elassieal probabihty current 
densities jQ{x,t) and Jc{x,t) are plotted for the values of m=l a.m.u, m=50 
a.m.u and m=720 a.m.u, corresponding to ctq = 10~ cm, u = 10'' cm/sec, C=50 
and X = I cm. The bold and the dashed curves represent the quantum and 
classical cases respectively. 
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Figure 2. The time variations of the quantum and classical probability current 
densities jQ(x,t) and Jc{x,t) are plotted for the values of C=100, C=10 and 
C=l by taking the mass= 1 a.m.u corresponding to (Tq = 10~^ cm, u = 10"* 
cm/sec, and X=l cm. The bold and the dashed curves represent the quantum 
and classical cases respectively. 



masses to be m = 1 a.m.u(H atom), m ~ 50 a.m.u and m = 720 a.m.u(C6omolecule). 
It is seen from Figure 1 that while for m = 50 a.m.u the disagreement between the 
quantum and the classical plots diminishes as compared to that for m = 1 a.m.u, a 
complete agreement is ensured from the masses around m = 720 a.m.u (Cgo molecule). 
In order to complement this line of study, another comparison is made between the 
plots of the quantum and the classical probability currents for a fixed mass, say, 
m = 1 a.m.u by varying the values of C ranging over C = 100, C = 10 and C ~ 1. 
Interestingly, it is seen from Figure 2 that while for C ~ 100, the quantum and 
the classical curves appreciably differ, this difference gradually diminishes with the 
decreasing values of C(i.e., as the departure from the minimum-uncertainty-product 
Gaussian wave function gets minimised) , with the difference becoming negligibly small 
as the value C = 1 is reached. 



Mass{a.m.u) 


{tc) (ms) 


{tg) (ms) 


{At)c (ms) 


{At)Q (ms) 


1 


68.684 


141.671 


16.592 


22.543 


5 


13.008 


14.575 


8.069 


9.523 


25 


12.881 


13.172 


7.940 


8.199 


50 


12.877 


13.025 


7.936 


8.064 
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100 


12.876 


12.947 


7.935 


7.999 


500 


12.876 


12.890 


7.935 


7.947 


1000 


12.876 


12.883 


7.935 


7.941 


5000 


12.876 


12.877 


7.935 


7.936 


10000 


12.876 


12.876 


7.935 


7.935 



Table 1. The comparisons between quantum and classical mean arrival times and their respective 
fluctuations are given for the different values of mass, corresponding to a fixed value of C = 1000, 
and the other relevant parameters being u = lOcm/s, X = 0.1cm, ctq = 10~*cm. 

Now, we come to a crucial aspect of this quantitative study; i.e., the probing of 
the range of masses over which an agreement emerges between the quantum and the 
classical mean arrival times, as well as between their respective fluctuations. For this 
we proceed as follows. We take three different values of C, viz. C = 1000, C = 500 
and C = 100, and, for any such given value of C, we vary the masses ranging from 
1 a.m.u(H atom) to the heavier molecules, say, biomolcculcs with molecular weights 
around lO'^ — lO^a.m.M (i.e., biomolcculcs comprising approximately 10-300 base pairs 
of DNA molecules, where 1 base pair « 650 a.m.u). 



Mass{a.m.u) 


{tc) (ms) 


{tg) (ms) 


{At)c (ms) 


{At)Q (ms) 


1 


115.021 


150.207 


22.686 


24.223 


5 


10.397 


11.324 


4.808 


5.654 


25 


10.281 


10.456 


4.711 


4.865 


50 


10.277 


10.365 


4.708 


4.784 


100 


10.276 


10.321 


4.707 


4.745 


500 


10.276 


10.285 


4.707 


4.715 


1000 


10.276 


10.281 


4.707 


4.711 


5000 


10.276 


10.276 


4.707 


4.707 



Table 2. The comparisons between quantum and classical mean arrival times and their respective 
fluctuations are given for the different values of mass, corresponding to a fixed value of C = 500, and 
the other relevant parameters being u = 10cm/ s, X = 0.1cm, ctq = 10"'* cm. 

Then, from the relevant computational results as given in Table 1, corresponding 
to C = 1000, it is seen that an appreciable difference between the quantum and the 
classical mean arrival times, as well as a significant difference between their respective 
fluctuations persist up to masses around 10"^ a.m.u, after which these differences 
gradually diminish. Eventually, these differences disappear beyond the mass range of 
10^ a.m.u(say, for the protein molecule such as cytochrome-c having the mass 12 x 10'^ 
a.m.u). It may also be noted that while the variations of both the quantities {t)c and 
{At)c as the mass changes saturate at the mass value of 10^ a.m.u, the corresponding 
variations of both the quantities (<)q and {At)Q with mass saturate around the mass 
value lO'* a.m.u. 



Mass{a.m.u) 


{tc) (ms) 


{tQ) (ms) 


{At)c (ms) 


{At)Q (ms) 


1 


201.172 


187.203 


28.556 


27.540 


5 


10.124 


10.321 


1.206 


1.697 


25 


10.037 


10.126 


1.064 


1.329 


50 


10.001 


10.002 


1.003 


1.128 


100 


10.000 


10.009 


1.000 


1.032 


500 


10.000 


10.002 


1.000 


1.006 
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Table 3. The comparisons between quantum and classical mean arrival times and their respec- 
tive fluctuations are given for the different values of mass, corresponding to a fixed value of C = 100, 
and the other relevant parameters being u = lOcm/s, X = 0.1cm, ctq = 10"'* cm. 

The results similar to that given in Table 1 are presented in Table 2 and Table 
3 for C = 500 and C = 100 respectively. While for C = 500, the agreement between 
{t)c and {t)Q, as well as between {At)cSLnd (Ai)Q emerge from the value of m=5000 
a.m.u(say, for the insulin molecule having m=5808 a.m.u), for C = 100, such an 
agreement is ensured from m=1000 a.m.u(say, for the peptide hormone molecule 
oxytocin having m=1007 a.m.u). 

It is, therefore, seen from these numerical computations that for a given set 
of values of the parameters u, X, and uo, the greater the value of C signifying an 
increasing deviation from the minimum-uncertainty-product Gaussian wave packet, 
the larger is the value of mass from which the quantum-classical agreement occurs for 
both the mean arrival time and its fluctuation. In other words, for the type of example 
studied here, by increasing the values of the parameter C, one can extend the range of 
mass values covering the heavier molecules for which appreciable disagreements can be 
found between the quantum and the classical values of the mean arrival time and its 
fluctuation. Thus, as regards the possibility of the relevant empirical studies, it would 
be interesting to probe the experimental rcalizability of the Gaussian wave packets 
that are characterised by large values of the parameter C. 

5. Concluding remarks 

An application of the specific quantum-classical comparison scheme underlying this 
paper is currently under consideration by using such a scheme for analysing the 
classical limit of quantum time distributions which are appropriately defined in 
the context of the harmonic oscillator potential - this study intends to use both 
the minimum as well as the non-minimum-uncertainty-product initial wave packets, 
including the case of the initial minimum-uncertainty-product wave packet having a 
specific width that corresponds to what is known as Schrodingcr's coherent state. 
Besides, further investigations are required to be pursued along, say, the following 
directions: 

a) The treatment presented in this paper uses essentially a Gaussian wave packet. 
It should, therefore, be interesting to study in terms of a suitably constructed non- 
Gaussian wave packet, or by using a superposition of wave packets, the extent to which 
the calculated quantitative results arc dependent on the form of the wave packet being 
strictly Gaussian. Such a study, also using forms of potential other than the linear one, 
would be particularly helpful for examining the empirical feasibility of tests related to 
the type of example discussed in this paper. 

b) Since the initial phase space distribution function used in our classical 
calculation is not uniquely fixed even if the position and momentum distribution 
functions are specified for a given wave function, it would be instructive to compare 
the quantitative results of this paper with the corresponding results for the choices of 
the initial classical phase space distribution function other than the simplest possible 
choice we have used for the given position and momentum distributions. As a special 
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case of such a choice, for a given non-minimum-unccrtainty-product Gaussian wave 
function, one may adopt the prescription given by Wigner[5^ for fixing the initial 
phase space distribution function to be used in our calculation. Studies along this 
line, based on the specific quantum-classical comparison scheme that has been invoked 
in the present paper, should be useful in throwing light on the extent to which the 
delineation of the mass values over which the quantum-classical agreement emerges 
for the mean arrival time and its fluctuation is sensitive to the choice of the initial 
classical phase space distribution . 

c) The quantum arrival time distribution used in this paper has been calculated 
specifically in terms of the probability current density. Comprehensive investigations 
are required in order to compare the quantitative results given in this paper with 
those obtained from various other schemes [17-21, 25-28] that have been suggested for 
defining the quantum arrival time distribution. This line of study would have an added 
ramification as regards the important issue concerning the possibility of discriminating 
between the different quantum approaches that have been proposed to compute the 
arrival time distribution. 

d) We note that the time-of-flight image method [34] has been widely employed 
for inferring the temperature of a cloud of trapped atoms in the experiments 
which involve the laser cooling of atoms. In this context, a quantitative study of 
the way the semi-classically computed time-of-flight distributions match with the 
corresponding quantum distributions in the limits of large mass and high temperature 
is of considerable interest [35]. It should, therefore, be relevant to take a fresh look at 
this issue by invoking the specific scheme for quantum-classical comparison that has 
been used in this paper. 
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